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Abstract 

We augment the hidden Stueckelberg Z' model by a pair of Dirac fermions in the hidden sector, 
in which the Z' has a coupling strength comparable to weak scale coupling. We show that this 
hidden fermion-antifermion pair could be a milli-charged dark matter candidate with a viable relic 
density. Existing terrestrial and astrophysical searches on milli-charged particles do not place 
severe constraints on this hidden fermion. We calculate the flux of monochromatic photons coming 
from the Galactic center due to pair annihilation of these milli-charged particles and show that 
it is within reach of the next generation 7-ray experiments. The characteristic signature of this 
theoretical endeavor is that the Stueckelberg Z' boson has a large invisible width decaying into the 
hidden fermion-antifermion pair. We show that existing Drell-Yan data do not constrain this model 
yet. Various channels of singly production of this Z' boson at the LHC and ILC are explored. 
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I. INTRODUCTION 



The standard model (SM) of particle physics has been blessed with her elegant way of 
giving masses to the weak gauge bosons by the Higgs mechanism. However, the crucial 
ingredient of this mechanism, the Higgs boson, is still missing. In addition, a scalar Higgs 
boson mass is not stable under perturbative calculation. It will receive an enormous amount 
of radiative corrections to its mass such that a delicate cancellation between its bare mass 
and radiative corrections is needed so as to obtain a mass in the electroweak scale - this is 
the famous hierarchy problem. An alternative way to give mass to an abelian U(l) gauge 
boson is known as the Stueckelberg mechanism. Although it is very difficult to give masses 
to nonabelian gauge bosons without sacrificing renormalizability within the Stueckelberg 
approach, it is worthwhile to study the consequence of this mechanism as an extension to 
the SM with extra abelian U(l) factors. 

Recently, Kors and Nath |l| showed that the SM extended by a hidden sector described 
by a Stueckelberg U(l)x and the gauge field C M associated with it can pass all the existing 
constraints from electroweak data as well as direct search limits from the Tevatron. Through 
the combined Stueckelberg and Higgs mechanisms, the SM SUl(2) x U(l)y gauge fields _B M 
and mix with the hidden sector gauge field C M . After rotation from the interaction basis 
(Cfj,, W^) to the mass eigenbasis (Z' Z M , A^), one obtains a massless state identified 
to be the photon 7 and two massive eigenstates which are the SM-like Z boson and an 
additional Z' boson. As long as the mixing is small, the Z' boson only couples very weakly 
to the SM fermions, and so it can evade all the existing constraints on conventional Z' 
models. The allowed mass range for the Stueckelberg Z' can be anywhere from 200 GeV to 
a few TeV j^j]. Typically, the mass of the Stueckelberg Z' is above the SM Z boson mass. 

In this work, a pair of hidden Dirac fermions is introduced in the Stueckelberg Z' model. 
Such a possibility has been mentioned in Ref. but its phenomenology was not explored. 
There could be various types or generations of fermions in the hidden sector, just like our 
visible world. Since the abelian U(V)x is assumed to be the only gauge group in the hidden 
sector and there is no connector sector between our visible world and the hidden one in this 
class of models, all hidden fermions in this sector are stable. 1 Thus the hidden fermion- 



1 This is in analogous to the pure QED case, muon does not decay into an electron plus a photon. 
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antifermion pair that we add in the model can be viewed as the lightest ones in the hidden 
sector, should there be more than one type of them. The SM fermions are neutral under this 
hidden U(l)x- Since this hidden fermion pair is stable, it can be the dark matter candidate 
of our Universe. 

In the next Section, we will present some details of the Stueckelberg Z' extension of the 
SM with an additional pair of Dirac fermion-antifermion in the hidden sector. In Section 
III, we discuss milli-charged dark matter. Treating the hidden fermion as our candidate 
of dark matter, we calculate its relic density and explore the parameter space allowed by 
the WMAP measurement. We also calculate the monochromatic photon flux coming from 
the Galactic center due to pair annihilation of these hidden fermions. In Section IV, we 
explore some novel collider phenomenology of the Stueckelberg Z' with the presence of the 
hidden fermion. Since the width of the Stueckelberg Z' is no longer narrow, compared to 
the scenario studied in {2], its phenomenology is rather different. Comments and conclusions 
are given in Section V. 



II. THE MODEL 

The Stueckelberg extension [l| of the SM (StSM) is obtained by adding a hidden sector 
associated with an extra U(l)x interaction, under which the SM particles are neutral. 2 We 
explicitly specify the content of the hidden sector: a gauge boson and a pair of fermion 
and antifermion x an d X- 

The Lagrangian describing the system is £stSM = £sm ^ st ' wnere 

£sm = - \w a ^ v W a » v - + D»$ - V(& $) 

+ ijrfD^f , (1) 

£st = ~\c, u C" + ixfDfx + \ (<V + M iC, + M 2 B,f , (2) 

D, = d^ + ig^WZ + igy^B^ (3) 

D* = d^ + igxQxC^, (4) 

where W^ u (a = 1, 2, 3), B^ u , and C^ u are the field strength tensors of the gauge fields W£, 



2 It was shown in Ref. [lj that the SM fermions are neutral under the extra U(l)x has the advantage of 
maintaining the neutron charge to be zero. 
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Bp, and C^, respectively. The SM fermions / were explicitly forbidden from carrying the 
U(l)x charges, as implied by Eq. ([3]), while the hidden fermion pair only carries the U(l)x 
charge, as implied by Eq. (J3J). One can show that the scalar field a decouples from the 
theory after gauge fixing terms are added upon quantization. 

The mass term for V = (C M , B^, Wjj) , after electroweak symmetry breaking ($) = 
v/y/2, is given by [l| 



/ 



l -V T MV = -~ (C M , B„ Wl 



Ml 



M X M 2 







M X M 2 M 2 + \g 2 Y v 2 -\g 2 g Y v 2 



-\929yv 2 



( a, \ 



B, 



\ -\g2gw 2 \g 2 2 v 2 ) \W*J 
A similarity transformation can bring the mass matrix M into a diagonal form 

fcA (z>\ 

B, =0 



(5) 



A 



T MO = diag(m|,, m|, 0) 



(6) 



The m 2 z , and m 2 z are given by 



m z> z 



Ml + M 2 + -(g 2 Y + g\)v 2 



±J(M? + Mf + \g 2 v 2 + ±g 2 v 2 ) 2 - (M 2 (g 2 + g 2 )v 2 + g 2 M 2 v 2 ) 



(7) 



The orthogonal matrix O is parameterized as 



O 



\ 



(8) 



-c e Sip cgc^ s e ) 

where = sin0, = cos</> and similarly for the angles ip and 9. The angles are related to 
the original parameters in the Lagrangian £stSM by 



X - + A M2 f ^ COS ' 

o = tan0 = — — , tan# = 

Mi g 2 



tan^ 



tanfl tan0m^ 



cos6[m 2 z , — mlyil + tan 2 9)] 
where m w = g 2 v/2. The Stueckelberg Z' decouples from the SM when — > 0, since 

M 2 



(9) 



tan i 



Mi 







tan ^ — > and tan 6* — > tan 6* v 



We note that the middle column is chosen to be different from that of Ref. [l| by an overall minus sign. 



where 9 W is the Weinberg angle. 

The interactions of fermions with the neutral gauge bosons before rotating to the mass 
eigenbasis are given by 

- = 92WI $ n ^f + 9yB, $rf^t + 9xC, x^Qxx , (10) 

where / denotes the SM fermions. The neutral gauge fields are rotated into the mass 
eigenbasis using Eq. (|S]), and the above neutral current interaction becomes 

- = $ f r \U l ,Pl + 4> Pr) Z' + U L P L + efrPn) Z, + eQ em A J ^ 



e*A M + e|Z M + e§,Z'l x , (11) 



where 



= gxQxi-ceSj,) , 
e| = 9xQ x ( s ^ c <P + ses^) , 

i 2 ^ 2 2 f-^y^ + g£) + s^s^- , 

92+9y4 V 2 27 2 

S " (12) 



fl-l + 9y c I 

In the above, we have used the relations 

r 3 Y 

e = 92 se = gyCtpCe and Q em = — + — , 

where Q cm is the electric charge operator. From Eqs. (fTTT) - (fl2|) . it is clear that in this class of 
model, the SM fermions interact with the hidden world through Z' and the hidden fermion 
interacts with our visible world through 7 and Z. In our computation, we assume the 
following input parameters at the electroweak scale [s| 

cwOz) = — 3 — , G F = 1.16637 x l(T 5 GeV~ 2 , m z = 91.1876 GeV , sin 2 W = 0.231 , 
128.91 

and the following three inputs related to the hidden sector 

5 = tan</>, gx, and Mz> ■ 

Since Q\ always enters in the product form gxQxi one can se t Qx to De unity without loss 
of generality. We derive from a cm , Gf, mz, and sin 2 9 W the values of 

-1/2 



y/4ira em , v = [V2Gf) , m w = m z yl — sin 6 W , and g 2 = 2m w /v 



We then fix the value of gy by the following equation 

e = = . 

\jgl + 9y4 

The other two angles 9 and ip are determined from the last two formulas in Eq. . 

It is clear from Eqs. f|TTl) - f|T2l) that the chiral couplings of the SM Z boson are affected 
by the mixing. In fact, even the mass of the Z boson is modified in this model, as shown by 
Eq. ([7j). It has been shown in Ref. that in order to keep the Z boson mass within the 
experimental uncertainty, the mixing angle must satisfies 



5 < 0.061 y/1 - {m z /M x f . (13) 

When 5 is small and mz> is large, M\ m m z > + 0(g2v). The constraint coming from the 
electroweak precision data is more or less the same as in Eq. flT3|) 

The limits obtained in Ref. also included the analysis from direct Z' production at 
the Tevatron. They showed that with the current Drell-Yan data, 

m Z ' > 250 GeV for 5 « 0.035 , 

m z > > 375 GeV for 5 « 0.06 . (14) 

If including the presence of a hidden fermion that the Stueckelberg Z' can couple to, the 
limit from direct Z' direction would be relaxed because of the smaller production rate into 
visible lepton pairs [2J. In Sec. IV, we will show that with a hidden fermion x fulfilling the 
dark matter constraint, the Z' would dominantly decay into the hidden sector fermion pair 
provided that mz> > 2m x . It would therefore entirely remove the constraint in Eq. (fT4"j) 
from direct production. On the other hand, if m z > < 2m x the Z' boson cannot decay into 
the hidden sector fermions, and so the constraint in Eq. (fill) stands. 

In the following numerical works, we will apply the constraints on 8 and m z > given by 
Eqs. (1131) and (fl4|) respectively, but when mz> > 2m x the latter constraint will be dropped. 

III. DARK MATTER PHENOMENOLOGY 
A. Milli-charged dark matter 

Milli-charged dark matter was first discussed by Goldberg and Hall 4|, motivated by 
the work of Holdom jlj] in which milli-charged particles in the hidden sector can interact 
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with the visible photon due to kinetic mixing between the visible photon and the hidden or 
shadow photon. Numerous constraints for the milli-charged particles, including accelerator 
experiments, invisible decay in ortho-positronium, SLAC milli-charged particle search, Lamb 
shift, big-bang nucleosynthesis, dark matter search, search of fractional charged particles in 
cosmic rays, and other astrophysical reactions, were summarized in jg] (see Fig. 1 of the 
first reference in jg]). Study of the constraints on milli-charged particles from neutron stars 
and CMB measurements were discussed in Refs. and 3 respectively. In summary, milli- 
charged particles of mass from MeV to TeV with a fractional electric charge (10~ 6 — lO^ 1 ) of 
a unit charge are still allowed. We note that integral charged dark matter was contemplated 



in 



9| and composite dark matter was studied in [lOj. More recently, PVLAS reported 



a positive signal of vacuum magnetic dichroism. It has been suggested [12] that photon- 
initiated pair production of milli-charged fermions with a mass range between 0.1 to a 
few eV and a milli-charge O(10~ 6 ) of a unit charge can explain the signal. However, this 
signal has not been confirmed by other experiments like the Q & A experiment For 
detailed analysis of various experiments of vacuum magnetic dichroism, we refer our readers 
to Ref.Q. 



B. Relic density and WMAP measurement 

The first set of processes we consider in our relic density calculation are 

XX -> Z'tZ^-t ff 

where / is any SM fermion. The amplitude for the annihilation x{Pi) x{P2) — * /(<7i) /fe) 
can be written as 



M = v(p 2 ) 7 M «(Pi) x M (<?i) 7 M [£lPl + £rPr) v(q 2 ) 



(15) 



where P L>R = (1 =F 7s)/2, and 



fXpOf ,xJl,r x Jl,r 

C 7 C -Vem t Z t Z , t Z' t Z' 



+ 2 + 

s — m z s — m Z i 



The differential cross section is given by 
da N f P f 



dz 327r s/3 A 



(16) 



Hi + m< + C + 2m^(s - 2m))) + Ami £ L £ R (s + 2m 2 ) 



(17) 
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where (3f tX = (1 — 4m| x / s) 1//2 , Nj = 3 (1) for / being a quark (lepton), t m = t — m 2 — m 2 = 
— s(l — (3f/3 x z)/2, u m = u — m^ — m 2 = —s(l+[3ff3 x z)/2, s is the square of the center-of-mass 
energy, and z = cos with the scattering angle. 

We also consider pair annihilation of XX into two neutral gauge bosons, 

XX -> W 2 with Vi,2 = 7, Z' (18) 
in our relic density calculation 4 . The differential cross section is given by 



|-2(2m 2 +m 2 yi ) (2m 2 +m 2 , 2 ) 



x "x^ 

2 i 2 , 2 N 



^X / \ ^x 



4 



u x t x 



(2 m 2 + + m^) (2 m 2 - - m^J ^9(2m x - m Vl - m V2 ) 

(19) 



where (3viV 2 — ^ 2 (l> m vJ s i m y 2 / s ) with A(a, 6, c) = a 2 + b 2 + c 2 — 2(ab + be + ca) is the 
Mandelstam function, t x and u x are given by t x = t — m 2 and u x = u — m 2 respectively, and 
S* is the statistical factor. We note that processes XX ~^ Hi ZZ are doubly suppressed by 
the small mixing angles and XX ~^ Z'Z' are either suppressed or forbidden by phase space, 
and therefore their contributions are negligible in the annihilation rates. 

The quantity that is relevant in the relic density calculation is the thermal averaged 
cross section (cry), where v is the relative velocity of two annihilating particles. In the non- 
relativistic approximation, v ~ 2 j3 x . To estimate the relic density of a weakly- interacting 
massive particle, we use the following formula [ljj] 

0.1 pb 



With the most recent WMAP 



n x h 2 _ 

(av) 



16] result on dark matter density 



Q CBM h 2 = 0.105 ±0.009 



where we have used the WMAP-data-only fit and taken f2cDM = ^matter — ^baryon, one can 
translate this WMAP data into 

(av) ~ 0.95 ± 0.08 pb . (20) 



We have ignored the channel \X ~ > J? Z' — > W + W which may contribute to certain extent. 
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FIG. 1: (a)-(b) are contours of av = 0.95 ± 0.16 pb (2 a range) in the plane of (gx, 8) for various 
mz> an d n^x- P &r t ( c ) shows the annihilation rate av versus m x with m/< = 500 GeV, gx = <?2> 
and 5 = 0.03. Part (d) shows the contour of av = 0.95 ± 0.16 pb (2 a range) in the (m x ,mz') 
plane. 

In estimating the annihilation rate during the freeze-out, we assume that the species has a 
velocity-squared v 2 ~ 0.1. To get a crude estimation, we ignore the thermal average and 
evaluate av directly. 

In FigOJa) and (b), we show the contours of av = 0.95 ±0.16 pb (2 a range) in the plane 
of (gx, S) for various input values of m z > and m x . We have included XX ~^ 1%', ZZ', and 
//, with / = u e , Vp, u T , e~ , \f ', r~, u, d, s, c, b, and t that are kinematically allowed. From 
Fig. 1(a) for m x = 60 GeV and mz> = 300 GeV, we can see that 6 ~ 0.03 and gx — 0.6 can 



give the correct amount of dark matter. Similarly, from Fig. 1(b) with the same m y 



60 



GeV and a larger mz> = 600 GeV, 5 ~ 0.03 and a slightly larger gx — 0.7 can also do the 
job. For comparison, we note that e ~ 0.31 and gi ~ 0.65 in the SM. Thus, the value of the 
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hidden Ux{1) coupling gx that we deduced from the WMAP measurement has the same 
order of the weak coupling g% of the visible sector. In FigQJc), we show the annihilation rate 
(tv versus m x for 5 = 0.03, gx = g-i, an d a fixed m Z i = 500 GeV. Clear resonance structures 
of Z and Z' are seen. In Fig{T](d), we show the contours of av in the (m x , mz>) plane. There 
are two branches: (i) the upper branch where m x < mz>/2 and the band relating m x and 
mz' is relatively wide; (ii) the lower branch where 2m x > mz> and the band relating m x and 
mz> is quite narrow. A narrow band implies the need of a fine-tuned relation between m x 
and mz' in order to give the correct dark matter density. In the latter branch, the Tevatron 
bound of mz> > 250 GeV for 5 ~ 0.03 has to be imposed. Therefore, the case of m x < m z >/2 
is more preferred theoretically. 

The hidden fermion x couples to the photon via the mixing angles cgs^, the value of 
which is about 0.9 x 0.03 ~ 0.03. Therefore, effectively the fermion \ "acquires" an electric 
charge of gxQx c s s <t>l e ~ 0.06, when its coupling to the photon is considered. Therefore, 
the range of m x ~ O(100) GeV and the size of effective electric charge ~ 0.06 implied by 
dark matter density requirement in our calculation are consistent with the constraints on 
milli-charged particles [6J]. 



C. Indirect detection 



If milli-charged hidden fermions like x an d X are the dark matter, their pair annihilation 
into 77, jZ, and jZ' in regions of high dark matter density, e.g. the Galactic center, can 
give rise to monochromatic 7-ray line that can reach our Earth for their indirect detection. 
The cross sections for these processes can be obtained from Eq. (fT9l) readily. The observed 
7-ray flux along the line-of-sight between the Earth and the Galactic center is given by [ljj] 

E) = crv^-^ f dsp 2 (r(s, ^)) , (21) 

d£/ 7 Anm x June of sight 

where the coordinate s runs along the line of sight in a direction making an angle ift from the 
direction of the Galactic center, dN 1 /dE 1 is the energy spectrum of the 7-rays, and v ~ 2/3 x 
is the relative velocity of the dark matter x an d X-, and the present value of v ~ 10 -3 . The 
flux from a solid angle AQ is often expressed as 

$JAQ,E) «5.6 x 10^ 12 ^ ( ^) (^^) 2 J(An)Ancm- 2 s- 1 , (22) 

dE 1 ypb j \ m x J 
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FIG. 2: The resulting photon flux from annihilation processes XX ~~ ^ 77; 7-^? an d jZ'. We have 
used typical values of J = 100, = 10~ 3 , and the present value of v as 10 -3 . 



with the quantity J(V0 defined by 

m 



8.5 kpc V0.3GeV/cm 3 



line of sight 



dsp 2 (r(s, ip)) 



(23) 



For the process of XX ~ > 77? we would have a mono-energetic 7-ray line with dN~/dE~ 



25{E 1 



Such a line, if observed, would be a distinctive signal for dark matter anni- 



hilation. Similarly, processes XX ~~ * 1% an d XX ~^ 1%' wm have a photon energy spectrum 



as dN-y/dEy rs 



m x (l 



m| 2 //4m^)). The contributions from these processes to the 



photon flux are shown in Fig. [2] with J = 100 and AQ = 10 3 . In this plot, we have taken 
a moderate value for J = 100 (averaged over 10~ 3 sr at the Galactic center). From Table 

n - 

7 of Ref. [15(1, we know the value of J varies from 2.166 x 10 (Kra profile) to 1.352 x 10 3 
(NFW profile) and to 1.544 x 10 5 (Moore profile). There are also cold dark matter profiles 
with dense spikes 17J near the Galactic center due to the accretion by central black holes 
that can give rise significant enhancement to the quantity J. With a rather conservative 
choice of J = 100, the flux of the 7-rays from the process XX ~~ *■ 77 is quite small due 
to double suppression of (e*) 2 . The process XX ~ ¥ 1% contributes at a somewhat lower 
flux level. The process XX ~^ 1%' can a l so contribute to the monochromatic photon flux, 



11 



provided that 2m x > m^. Since this process is only suppressed by one power of the mixing 
angle, it could be more substantial than the doubly-suppressed process XX ~* 77- Note 
that since 2m x > m^, the Tevatron bound Eq. (ll4j) of mz> > 250 GeV for S ~ 0.035 must 
be enforced. When kinematics allowed, the photon flux from this process is three orders 



of magnitude higher than that from XX ~^ 77- F° r com 



3arison, we note that the photon 



flux from the neutralino pair annihilation x°X° ~~ * 77 1S| nas been estimated to be about 
1.5 x 10 -14 (2 x 10~ 13 ) cm _2 s _1 if the neutralino is a Higgsino-LSP (Wino-LSP), using the 



same moderate value of J = 100 



3, 



20]. 



The expected sensitivities for the new Atmospheric Cerenkov Telescope (ACT) experi- 



ments such as HESS |2l| and VERITAS |22| are at the level of (10" 14 - 10~ 13 ) cm" 2 s" 1 
with an angular coverage of about 10 -3 . They are sensitive to the 7-rays from a few hundred 
GeV to TeV. On the other hand, the Gamma-ray Large Area Space Telescope (GLAST) 



experiment [23j due for launch in Fall this year, can probe 7-rays from 20 MeV to 300 GeV, 
but at a lower sensitivity level about 2 x 10~ 9 — 10~ 10 cm~ 2 s _1 . From Fig. [21 there is 
a small range of m x ( m x < 100 GeV) such that XX 77 contributes at a level larger 
than 10 -14 cm" 2 s -1 . The process XX ~^ 1% contributes at a level below the sensitivities of 
all these experiments for all ranges of m x , whereas the process XX ~ ¥ 1%' can contribute 
at a much larger flux and it is indeed above the sensitivity level of ACT experiments for 
m x < 600 GeV. Since GLAST can only be sensitive to 7-rays of 300 GeV or less with lower 
sensitivity, it is hard to detect the 7-rays from the lighter milli-charged dark matter. For 
heavier milli-charged dark matter, the 7-rays can be above a few hundred GeV and thus 
above the sensitivity reaches of HESS and VERITAS. 

The continuum 7-ray background from astrophysical sources near Galactic Center is 
largely an unsettled issue due to astrophysical uncertainties. There have been data showing 
excess of 7-rays in different energy regimes near the Galactic center. The EGRET experiment 



24( has reported an excess of 7-rays in the region of the Galactic center, including the galactic 



longitude and latitude position at / = 0° and 6 = 0° degrees. The level of excess is above 
the expectation of primary cosmic rays interacting with interstellar medium. The EGRET 
excess region is around 1 GeV. However, there may be some other unknown sources of 7-rays 
around the Galactic center. It is hard to establish the fact that the excess is due to dark 
matter annihilation, because the excess does not have specific features. This is in contrast 
to the monochromatic 7-ray flux, which is a clean signature of the dark matter annihilation. 
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In the Galactic center region, excesses of 7-rays were also reported by VERITAS 25J in 



the range above 2.8 TeV and by CANGAROO collaborations [26] in the range of 250 GeV to 
1 TeV. Such excesses are also hard to be explained by conventional dark matter candidates. 
The HESS Collaboration also had a measurement of TeV gamma rays from the Galactic 
Center 27|, which is, to some extent, in disagreement with the CANGAROO results. It 



was pointed out [28j that this TeV 7— ray excess is likely to be of astrophysical origin and 
thus it constitutes a background for detecting dark matter annihilation. The origin of these 
backgrounds may be due to violent acceleration of cosmic protons and other particles by the 
chaotic magnetic fields near the Galactic center black hole 29j . After escaping the black hole 
environment and fly off into the interstellar medium, these extremely high energy protons 
collide with low energy protons (hydrogen gas) to form pions, which subsequently decay into 
high energy 7-rays that can radiate in all directions. 

Due to its unknown astrophysical origin, it is hard to establish accurately the true contin- 
uum 7-ray background which could be used for comparison with dark matter annihilation. 
Thus, using the continuum 7-ray signal is difficult to provide strong evidence for dark mat- 
ter, unless the dark matter annihilation rate is very large. On the other hand, provided that 
the annihilation cross section is sufficiently large, the monochromatic photon line would 
be a "smoking gun" signal for dark matter annihilation, since the energy of the 7-ray is 
uniquely determined by the mass of the milli-charged dark matter (and the Z' mass as 
well in the XX ~ * 1%' channel). Nevertheless, the EGRET and HESS continuum back- 
grounds still pose a serious challenge to detecting the monochromatic photon line due to 
dark matter annihilation in the Galactic center region 28j. It was shown in Ref. 28j that 
in order for a photon line to be detected above the continuum background, the quantity 
(av/\0~ 29 cm 3 s _1 ) J AQ must be larger than 10 — 100. This implies that the photon flux 
to be larger than 1.9 x (TeV/m x ) 2 x (10~ 14 — 10~ 13 ) cm -2 s -1 . From Fig. [21 it is easy to see 
that for m x between 150 and 300 GeV, the photon flux in the XX 1%' channel is within 
the detectability level. 

One may also give a rough estimate for continuum photon flux arises from the milli- 
charged dark matter annihilation into light quark pairs. The continuum photon spectrum 
mainly comes from the light quark fragmentation into neutral pions, which subsequently 
decay into secondary photons. The differential spectrum dN 1 / dE^ can be obtained by Monte 
Carlo event generators, e.g. PYTHIA, and then parameterized as a quark fragmentation 
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function. We can use Eq. (l2Tj) with dN^/dE^ given by a fragment at ion- like function 30[: 



dN 

—1 =r]x a exp(b + cx + dx 2 + ex 3 ) , (24) 
dx 

where x = E 7 /m x and for a light quark, e.g. u or d quark at an energy of 500 GeV, rj — 1, 
a = —1.5, b = 0.047, c = —8.7, d = 9.14, and e = —10.3. These coefficients depend only 



mildly on the energy of the light quarks 30j . Putting all these factors together, we estimate 
the integrated photon flux with E 7 > 1 GeV to be of the order of 10 -10 (10~ n ) cm~ 2 s" 1 for 
m x = 100 (500) GeV. It is at most around or slightly below the sensitivity level of GLAST. 
Since the VERITAS and HESS experiments are sensitive to higher energy and the above 
spectrum Eq.( j24|) falls off rapidly as x increases, their integrated photon fluxes are at least 
an order of magnitude smaller than that of GLAST. Despite challenging by the uncertain 
astrophysical backgrounds, this continuous secondary photon spectrum together with the 
monochromatic photon line from milli-charged dark matter annihilation can be probed by 
the next generation of 7-ray experiments. 

IV. COLLIDER PHENOMENOLOGY 

Phenomenology of the Stueckelberg Z' with the presence of the hidden fermion- 
antifermion x an d X that the Z' can decay into is quite different from the one studied 
before in Refs. [ll, |2j. 

The partial width of Z' into a SM fermion pair // is given by 

r<*W/) = ^M(# + #) (l-g) + .efe*g] (25, 
and into hidden fermion pair XX is simply 

r (z ^«)^ mz , £ |, 2 (i + 5). (2 6) 

Here, (3f >x = (1 — Amj x /m 2 z ,) 1 / 2 . The total width of Z' is evaluated by summing over all 
partial widths, including both the SM modes and the hidden mode. We show in Fig. [3] the 
various branching ratios for Z' as a function of its mass with the following inputs gx = 92, 
5 = 0.03, and m x = 60 GeV. Since the mixing angle is so small (5 = 0.03), the Z' is 
mainly composed of the C M boson of the hidden sector. Hence, the Z' dominantly decays 
into the hidden sector fermion pair while it has only a small fraction of 10 -3 into visible 
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FIG. 3: Branching ratios for Z' with gx = 92, $ = 0.03, and m x = 60 GeV. 



fermions. The strategy for the search of this Z' would be very different from all the previous 
conventional Z 1 models including the hidden Stueckelberg Z' studied in 

A3. 

Before we explore for the possible collider phenomenology of the Stueckelberg Z' boson 
and the hidden sector fermion x, we have to make sure that the new particles and the hidden 
sector interactions will not upset the existing data. 



A. Constraints from invisible decays of Z and quarkonia 

Firstly, the SM Z boson that is observed at LEP would decay into a pair of hidden 
fermions xXi giving rise to additional invisible width other than the neutrinos. Because of 
the mixings among the three neutral gauge bosons, the Z boson can couple to the xX P & i r y i a 
the mixing angle s$. We have calculated the partial width of Z — > xX f° r 9x — <?2, 8 — 0.03 
(consistent with the dark matter requirement), and m x = — 45 GeV. The partial width is 
about 0.24 MeV^which is much smaller than the uncertainty 1.5 MeV of the invisible width 
of the Z boson 3|. Even if we allow a larger mixing angle 5 = 0.061 (its maximum value 
allowed by Eq. (fTBl ). the invisible width of Z would be at most 1 MeV, which is still within 
the la uncertainty of the data. If the mass of x is beyond half of the Z boson mass, the 
invisible width of the Z boson would not constrain the model. 

The hidden fermion x can a l so couple to the photon via the mixing angles cgs^, the 
maximum of which is about 0.9x0.03 ~ 0.03. Therefore, effectively the fermion x "acquires" 
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an electric charge of gxQx c 8 s 4>/ e ~ 0-06 when its coupling to the photon is considered. If 
X is very light, of the order of MeV, it could be produced in J ftp and T decays as invisible 
particles. Constraints on invisible decays of J lib and T exist (for a comprehensive review 



on constraints on light dark matter: see Ref. 



3l||). A very recent update on the T(15') 



invisible width is given in Ref. [32j. The invisible widths of J ftp and T are respectively 

B(J/ip -> invisible) < 7 x 10~ 3 and B(T(1S) -> invisible) < 2.5 x 10~ 3 . 

However, the partial width of J/ip into XX is suppressed by the milli-charged factor of 
(0.06) 2 relative to the partial width into e~e + . Thus B(J/ip — > xx) ~ (0.06) 2 x B(J/ip —> 
e~e + ) ~ 10~ 4 , which is well below the above limit. The situation for T invisible decay 
is very similar: B(T(1S) — > yy) « (0.06) 2 x £?(Y(1S') — > e~e + ) ps 10~ 4 , which is also 



safe. Indeed, a recent study |33j using 400 fb" 1 luminosity collected at the T(4S'), the B- 
factory can limit B(T(1S) — > invisible) ^ 10~ 3 . There are also other decays modes, such as 
J/ip or T — * 7 + invisible, but it is straightforward to check that with an effective charge 
of 0.06 the experimental limits of these radiative invisible modes do not constrain the model 
at all. If the mass m x is above 5 GeV, it has no constraint at all from these invisible decays 
of the quarkonia. 



B. Constraint from singly production of Z' at LEPII 



Singly production of the Z' at LEPII is possible via e~e + — * jZ' followed by the invisible 
decay of the Z' . This process is very similar to the SM process e~e + —* ^Z — * ^vv. The 
differential cross section for e~e + — > jZ' is given by 

da(e~e + -> 1 Z') f3 z >e 2 Ql 



( e L 2 e R 2 



1 

ut 



r 



u 



2sm 2 z , 



(27) 



dcosQ 327rs 

where 6 is the scattering angle of the photon, t = —s0z>(l — cos6)/2, u = —s/3z>(l + 
cosG)/2, and f3z> = (1 — m z ,/s). We show the production cross section at LEPII energy 
y/s = 205 GeV in Fig. H]as a function of my. Since the Z' would decay into invisible yy, 
the signal of which would be a mono-photon. The recoil mass spectrum would then indicate 
the mass of the Z 1 . In the figure, we also show the 95% C.L. upper limits on mono-photon 
production as a function of the missing mass obtained by DELPHI 34j . A small mass range 
of Z' , 180 GeV < m Z i < 200 GeV, is disfavored by the data. However, one has to be cautious 
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FIG. 4: Comparison with the DELPHI data on the mono-photon production. The theory prediction 
is for gx = <?2 and 5 = 0.03. 

in this relatively soft photon region where large theoretical uncertainties are expected to be 
important. 



C. Drell-Yan production of Z' at the Tevatron 

The production cross section of Z' followed by the leptonic decay at the Tevatron is given 

by 

" q=u,d,s,c 

(28) 

where y^i = 1960 GeV, r = m|,/s, IV is the total width of Z' given in Eqs. (|25|) and 
(I26p . and e^' fl can be found in Eq. Q12p . This Drell-Yan cross section for the Z' boson 
is plotted in Fig. [5j where the 95% C.L. upper limits on cr(Z') ■ B(Z' — > e~e + ) from the 



CDF preliminary results [35(] are also shown. It is clear that the present CDF limits do 
not constrain the model at all. This is in sharp contrast to the results studied in Ref. 
because the Z' boson that we consider here has a very small branching fraction into charged 
lepton pair. The Z' boson would decay preferably into the hidden sector fermions instead 
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FIG. 5: Drell-Yan cross sections pp — > Z' — > e e + versus mz' for gx = <?2 an d 5 = 0.03. We also 



show the 95% C.L. upper limits on <t(Z') ■ B(Z' — > e e + ) from the CDF preliminary results 

of visible particles. On the other hand, the Stueckelberg Z' in Ref. only decays into the 
SM particles. In our case the Z' only has a branching ratio of ~ 10™ 4 — 10~ 3 into leptonic 
pairs, and it would not be easily detected in the Drell-Yan channel. Neither the hadronic 
decay modes of Z' can afford it to be detected. Even in the future runs of the Tevatron with 
a sensitivity reaching the level of 10 -3 — 10 -2 pb, it is still not possible to detect this kind 
of Z' boson through the Drell-Yan channel. 



D. Singly production of Z' at LHC and ILC 

Perhaps one has to rely on the invisible decay mode of the Stueckelberg Z' of this model 
to identify its presence. Here we calculate the predictions of singly Z' production at the 
LHC and ILC. Other than the Drell-Yan process that we have considered, the next relevant 
process to probe this invisible Z' is via qq — > gZ' followed by Z' — > \Xi which gives rise to 
monojet events. The subprocess cross section can be easily adapted from Eq. (1271) : 

dajqq -> gZ') = p Z '9 2 s / qL 2 qR 2\ ± 
dcos6* 72tts \ z ' z ' > u< 



i 2 + u 2 + 2sm%, 



(29) 
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FIG. 6: Production cross section for the process pp — > j + Z' followed by invisible decay of the Z' 
with gx = 92 and 5 = 0.03 at the LHC. We imposed pxj > 20 GeV and |r/j| < 2.5 on the jet. 

Other cross channels, e.g., qg — > qZ', can be obtained from Eq. (1291) using the crossing 
symmetry. The branching ratio BiZ' — > XX.) is very close to unity. We show in Fig. [6] 
the production rate of monojet events versus with gx = gi and 5 = 0.03 at the LHC 
under the jet cuts of pxj > 20 GeV and \r]j\ < 2.5. The qqZ' coupling is suppressed by the 
small mixing angle, the same as in the Drell-Yan process, but unlike the Drell-Yan process, 
this monojet amplitude is suppressed by only one power of the mixing angle instead of two. 
Therefore, the rate is not negligible. Also, the true SM background for monojet is rather 
rare. Thus, the monojet event actually signals the presence of such an invisible Z' . 

Another place to detect such an invisible Z' is at the ILC with the process e~e + — * ^yZ' — > 
IXXi which we have considered above for the mono-photon limits from LEP. We extend the 
energy to 0.5 — 1.5 TeV and calculate the event rates for the mono-photon final state. We 
show in Fig. [7|the production rates at y/s = 0.5, 1, 1.5 TeV with g x = #2 and 5 = 0.03. 
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FIG. 7: Production cross section for the process e~e + — > 7 + Z' followed by invisible decay of the 
Z' with g x = 92 and 6 = 0.03 at ILC (y/s = 0.5, 1, 1.5 TeV). We imposed £ 7 > 10 GeV and 
I cos 7 I < 0.95 on the photon. 

V. CONCLUSIONS 



We have proposed an extension of the Stueckelberg Z' standard model by adding a pair 
of fermion and antifermion in the hidden sector, which has only a U(l)x symmetry. The 
stability of the hidden fermion pair with its weak sized interaction makes it a suitable dark 
matter candidate with a correct amount of dark matter density. We have calculated the 
photon flux from the Galactic center due to the annihilation of this milli-charged dark 
matter. If 2m x < mz>, there is only a small range of m x that the photon flux is above the 
sensitivity level of the future 7-ray experiments. However, if 2m x > mz> there is a wide 
range of m x that the photon flux is above the sensitivity level. The collider phenomenology 
may be different from those studied in Ref. Q], because the dominant decay of the Z' 
is into the invisible XX if kinematically allowed. In this case, the present Drell-Yan data 
cannot constrain the model at all. We have proposed the monojet signal at the LHC and 
the mono-photon signal at the future ILC to probe this invisibly decaying Stueckelberg Z' 
boson. 
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We close with some comments. 



1. Since only a Ux(l) symmetry is assumed in the hidden sector, each hidden fermion is 
stable against decay. Therefore, if we assume more fermion pairs in the hidden sector, 
their relic densities are additive. Thus, a larger coupling constant is needed to ensure 
larger annihilation cross sections. One can also consider multiple hidden Stueckelberg 
U(l) extension of the SM. We refer to Ref.jjJ for the discussion for this possibility. 

2. When rnz' < ^^x^ the decays dominantly into visible particles. It can be easily 
detected in the Drell-Yan channel. The existing data constrains the model, as given 



by Eq. (1T4")) originally obtained by the authors in Ref.j^]. Photon flux from pair 
annihilation of XX ~^ 1%' a t the Galactic center is also within reach at the next 
generation of 7-ray experiments. 

3. When mz> > 2m x , the Z' decays dominantly into invisible XX- The present Drell-Yan 
data cannot constrain the model, neither can the invisible decays of J/ip and T for 
a very light X- However, the mono-photon production limits obtained by DELPHI 
disfavors a small range of 180 GeV < m z , < 200 GeV. We anticipate that in the 
future ILC the missing mass spectrum can efficiently constrain this type of invisibly 
decaying Z' boson. 

4. The hidden fermion appears to have a milli-charge as it acquires a small effective cou- 
pling to the photon through the mixing induced by the combined Higgs and Stueck- 
elberg mechanisms. With a mass of 0(100) GeV and an effective charge 0.06 of a 
unit charge, the hidden fermions are consistent with the existing constraints on milli- 



charged particles [6J. As milli-charged particle is of very recent interests [12j, an update 
on the terrestrial and astrophysical constraints on this hidden milli-charged particle is 
desirable. 
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Note added. Stueckelberg Z' extension with kinetic mixing has been studied recently in 



361 ] . Wherever overlaps in the parameter space, the authors in 36] found agreements with 



the analysis presented in our work. 
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